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The effect of external fluctuations on the formation of 
spatial patterns is analysed by means of a stochastic Swift- 
Hohenberg model with multiplicative space-correlated noise. 
Numerical simulations in two dimensions show a shift of the 
bifurcation point controlled by the intensity of the multiplica- 
tive noise. This shift takes place in the ordering direction (i.e. 
produces patterns), but its magnitude decreases with that of 
the noise correlation length. Analytical arguments are pre- 
sented to explain these facts. 



du (a?, t) 
dt 



= ru-[\7 2 + kl] 2 u-u 3 + 'n(x,t) (1.1) 



I. INTRODUCTION 

A large amount of experimental spatially-extended sys- 
tems exhibiting non-equilibrium transitions controlled by 
the environment exist (for a review, see (|). These tran- 
sitions usually occur in such a way that the system de- 
parts from an initial homogeneous state when an ex- 
ternal control parameter surpasses a certain threshold 
value. Frequently, this instability leads to the appear- 
ance of non-equilibrium spatiotemporal dissipative struc- 
tures, which subsist as long as the external driving stress 
which has produced the transition persists [Q. This is 
the case, for instance, in the hydrodynamic convective 
structures controlled by the temperature difference be- 
tween the plates of a Rayleigh-Benard cell |j||l| , or in the 
appearance of transverse structures in a laser beam con- 
trolled by the pumping rate [HQ. Usually, most of the 
features present in these pattern-forming processes are 
satisfactorily captured by a model equation which was 
introduced by Swift and Hohenberg Jt| with the initial 
objective of describing the effect of fluctuations at the on- 
set of convection in the Rayleigh-Benard cell mentioned 
above. Since then, the Swift-Hohcnberg (SH) equation 
(and proper modifications of it |^|| ) has proved its great 
usefulness in this field, not only in a hydrodynamical con- 
text [fL0l|ll|]. The equation reads: 



Here r is the bifurcation parameter (the temperature gra- 
dient in the Rayleigh-Benard case, or the pumping rate 
in the laser case), which controls the appearance of a pat- 
tern of characteristic length of order k^ 1 . The field r\ is a 
random function of space and time, representing internal 
noise of the system. Usually this term is statistically de- 
scribed by a gaussian probability distribution with zero 
mean and correlation 

(r)(x,t) r)(x',t')) = 2e5(x-x') 5(t-t') (1.2) 

This white noise accounts for hydrodynamic thermal fluc- 
tuations in the convective instability, or for spontaneous 
emission in the optical case. Its (dimensionless) strength 
e is very small, so it has no qualitative effect on the pitch- 
fork bifurcation exhibited by this model. 

External fluctuations, on the other hand, are known 
to have non-trivial qualitative effects on the behaviour 
of non-linear dynamical systems. In particular, zero- 
dimensional systems (i.e. systems with no spatial de- 
pendence) have been known for more than a decade to 
exhibit transitions induced by external noise (see 
for an extensive review). In the last few years, inter- 
est on studying the influence of external fluctuations on 
spatially-extended systems has grown [^3|-p7[. This in- 
fluence is likely to be more relevant in these systems than 
in the previous homogeneous ones, due to the existence 
of symmetry-breaking effects in the transitions occurring 
in them. 

In the particular case of the SH model, Elder et al. 
ijlij observed a disordering effect of the additive noise 
present in Eq. (1.1), when it is no longer considered 
to be internal, so that its intensity e can be arbitrarily 
large. Hence, the stripe structure (roll structure, in the 
hydrodynamical terminology) suffered a transition from 
a smectic (large regions of parallel stripes) to an isotropic 
(disordered system with short-range order) regime as the 
intensity of the additive external noise increased. 

Yet another possible (even more reasonable) source of 
external noise exists, namely the control parameter of 
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the bifurcation. The fact that this parameter is directly 
related to an external constraint imposed by the observer 
makes it likely to be affected by fluctuations. Let us now 
denote by r the mean value of this fluctuating control 
parameter. Now the SH equation reads: 



du (x, t) 
dt 



[r + i{x,t)\ u- [V 2 



t) (x, t) 
(1.3) 



Hence the random field £ is a zero-mean external 
multiplicative noise (chosen gaussian) whose correlation 
should not in principle be assumed to be white. We shall 
choose here a noise colored in space and white in time: 



(£(*,<) £(£',*')> = 2D 



A 



5(t-t') (1.4) 



A is the correlation length of the noise. This model was 
studied in Ref. pl| in the limit case A — ► (noise white 
also in space) and nontrivial effects were indeed found 
on the original bifurcation, which was shifted into the 
pattern region as the intensity of the external noise in- 
creased. Hence the role of this external noise is, at least 
for some values of its intensity, opposite to that of the 
additive case. A linear stability analysis of the first and 
second moments of the relevant variable of the system 
leads to analytical results which are in agreement with 
simulations pq , p!q| . Now we are interested in the influ- 
ence of the correlation length of the noise (which in some 
other problems is known to reduce the effective value 
of the noise intensity |^8|,^9|). The mathematical treat- 
ment of the Fokker-Planck equation to study spatially- 
extended systems with multiplicative colored noise is also 
presented. 

This paper is organised as follows. Next Section 
presents the results obtained by numerical simulation of 
the model. A linear stability analysis of the structure 
function of the system is made in Section [II, where a 



comparison with the previous numerical results shows 
good agreement. Section IV contains a generalisation 



of the linear stability analysis treatment to higher-order 
statistical moments. Finally some conclusions are stated. 
An Appendix contains some details on the Fokker-Planck 
approach to the problem. 



II. NUMERICAL ANALYSIS OF THE FULL 
NONLINEAR MODEL 

The non-equilibrium, nonlinear stochastic problem 
presented in the previous Section does not admit an exact 



analytical study. Hence our first approach to the prob- 
lem is a numerical one, and constitutes the best way of 
obtaining a first insight into the effects of external fluc- 
tuations in the pattern-forming instability developed by 
the SH model. 

The behaviour of the model in the absence of multi- 
plicative noise is well known. When e is small, the ho- 
mogeneous state u = is stable for a negative value of 
the external control parameter r. For r ~ 0, a bifur- 
cation takes place from this homogeneous situation to 
an inhomogeneous state composed of stripes. The in- 
fluence of an uncorrelated fluctuating control parameter 
in this non-equilibrium spatially-extended transition was 
preliminarly analysed in Ref. |15fl , revealing a nontrivial 
ordering effect. Here we extend the numerical analysis to 
a colored case. 



A. Algorithm 

We begin by discretizing space in a regular two- 
dimensional square lattice with L x L cells of size Ax. 
Now the SH model can be written in the following gen- 
eral form: 

^=/<(«)+fltf («)&(*)+*(*) (2-1) 

with 

< Vi (t) 7fc(f) > = 2 Sii 6(t - (2.2) 



<&(*)&(*') > = 2 Di-j 5(t-t') . 



(2.3) 



Cells are named with one index and repeated indexes are 
summed up. The deterministic force / and the coupling 
function g are in our particular case: 

fi (u) = r Ul - (V 4 u) i - 2 (V 2 u) i - m - u? (2.4) 



g lj (u) = Ui 5 tj (2.5) 
The discrete Laplacian operator is defined as 
1 



(Ax) 



2 Z-^i 



(u n - Ui) 



(2.6) 



where the sum extends over the set of nearest neighbours 
of site i. Now dynamical evolution is discretized in time 
(let At be the time integration step) and an algorithm 
can be developed from standard techniques |po|| . Up to 
O (At 3 / 2 ), it reads 
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Ui {t + At) = Ui {t) + fi (u(t)) At + gij (u(t)) X 3 

+ Y * + \lB 9ki X, X l + l^i X 3 Y k (2.7) 



2 du k 



2 du k 



where Yj is a gaussianly distributed random number with 
zero mean and variance 2Ate/Ax d to be placed at site 
j. Xj = \J2 At Zj, where Zj is a random field correlated 
in space, whose gen erat ion procedure will be described 
later. According to ( |2.5| ), the algorithm becomes in our 
particular case 

Ui (t + At) = u t (t) + fi (u{t)) At + u^X, 

+ ^ Ul (t)Xf +Y i + ~X i Y i (2.8) 

(now repeated indexes are not summed up). The last 
term in this expression can be replaced by its statistical 
average (which is zero here), since it is known that in 
this case the resulting algorithm represents a stochastic 
process having the same statistical properties . Hence 
the algorithm which we use is finally: 



Ui (t + At) = Ui (t) + fi (u{t)) At + u t (t)Xi 

1 



+ -Ui{t)Xf + Yi 



(2.9) 



The space-correlated random field Z will be generated 
from the following relation: 



Z(x) 



exp 



— A 2 V 2 
2?r 



W(x) , 



(2.10) 



which is written in continuum space. W (x) is a gaussian 
white random field of intensity D: 



(W(x) W(x')) = DS(x-x') 



(2.11) 



Z (x) can be easily generated in Fourier space, w here W 
becomes an anti-correlated field . Definition ( [2.10] ) is 
chosen in such a way that random field Z (x) (and hence 
multiplicative noise £ (if, t)) have a well-defined correla- 
tion length A, as can be seen by computing the space 

correlation function D (^ x ~^ ^ , which can be done in a 



straightforward way by writing ( 2. 1C ) in Fourier space. 
The result is: 
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(2.12) 



B. Analysis of the transition 

In order to anal yse the p attern-forming bifurcation ex- 
hibited by model (2.1-2^), we shall define the following 
quantity: 



J(t) 



u 2 (x, t) dx 



(2.13) 



where the statistical average is made over realizations of 
both the internal and external noises. In the Rayleigh- 
Benard case, the density of this quantity (j(t) = J(t)/V) 
corresponds to the density of heat flux due to convec- 
tion from the lower towards the upper plate of the cell. 
Hence its value is zero in the homogeneous state (no con- 
vective rolls) and non-zero in the structured convective 
phase, increasing linearly with the control parameter r, 
as shown clearly by experiments. This behaviour is re- 
covered by numerical simulations of the model. In the 
absence of external noise, stripes appear for values of r 
greater than r ~ 0, as discussed earlier. This can be seen 
in Fig. |l|, where the steady-state value of the heat flux is 
plotted against the control parameter of the system. It 
should be noted that the existence of a small but non-zero 
internal additive noise, along with the fact that simula- 
tions are performed on small finite systems, produces a 
rounding-off of the otherwise sharp transition described 
above. Periodic boundary conditions are considered, and 
Ax = 1 is chosen. Given this lattice spacing, a value 
of At = 0.01 happens to be enough to ensure stability 
of the algorithm. The intensity of the additive noise is 
taken to be equal to 10~ 3 in all cases. On the other hand, 
the value of the wavenumber fco is chosen in each case so 
that each convective roll is described by 10 lattice cells. 
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FIG. 1. Steady convective flux vs. control parameter r in 
the absence of external noise (circles), with a white multi- 
plicative noise of intensity D — 0.1 (squares) and with a mul- 
tiplicative noise of intensity D — 0.1 and correlation length 
A == 0.71 (triangles and stars). System size is 40 x 40 except 
for the stars, which correspond to 30 x 30. 

In the presence of a non-zero external multiplicative 
uncorrelated (A = 0) noise the bifurcation is shifted to 
the conducting (r < 0) region, which leads to the ap- 
pearance of roll patterns in otherwise subcritical regions 
(squares in Fig. l]). The amount of the shift is decreased 
back by the introduction of space correlation (A ^ 0) in 
the external fluctuations, as shown by the empty-triangle 
curve in Fig. El In this last case the simulation is also 
performed for a smaller value of the system size, but no 
finite-size effects are encountered (asterisks in the same 
figure) . 
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FIG. 2. Relative fluctuations of the steady convective heat 
flux vs. control parameter r. Values of the parameters are the 
same as those in Fig. 1. 

In analogy to studies of equilibrium phase transitions, 
the bifurcation shift found above can also be observed 
from the computation of the relative fluctuations of the 
steady heat flux, which are expected to exhibit some sort 
of singular behaviour in the transition point. Indeed, if 
one computes the relative fluctuations of J as 



of neither J nor its relative fluctuations x gives a hint 
on the effect of multiplicative noise on the nature of the 
bifurcation (i.e. on the order of the transition). Further 
analytical work would be needed in order to clarify this 
point. 

In conclusion, simulations show that stripe patterns 
can be favored by multiplicative noise, this effect being 
diminished by correlation of the noise. Figure || shows 
three different situations corresponding to the states 
marked A, B and C in Fig. [j]. Cases B and C are patterns 
for r < and hence favored by multiplicative noise. 




(b) 




V 



(j 2 )-(jy 



(2.14) 



then the transition point is characterised by a maximum 
value of this quantity. This is observed in Fig. where 
the transition shifts, which were already present in Fig. 
|l|, are clearly observed. Again no finite-size effects arc 
found. On the other hand, unfortunately the analysis 
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FIG. 3. Patterns corresponding to the three situations 
named A, B and C in Fig. 1. 



C. Structure function in the convecting phase 

In order to analyse more deeply the lack of finite-size 
effects in the convective heat flux above threshold, we 
have computed the stationary spherically-averaged struc- 
ture function in a noise-favored convecting state (with a 
negative r but a supercritical value of D). This func- 
tion is defined as the Fourier transform of the correlation 
function of the system, 

G(x,t) = — / dx' (u(x',t) u(x' + x,t)) (2.15) 

v JV 

Hence, the structure function indicates any periodicity 
in the system, so that it is a good way of characterising 
a pattern. Let us consider a finite system of volume V. 
The Fourier transform u^t) of the field u {x, t) is defined 
by: 



t 1 



(2.16) 
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FIG. 4. Structure function for the same convecting state 
with different system sizes. 

Simulations show that this function behaves in the 
same way for the whole parameter range, presenting a 
unique maximum related to the periodicity of the pattern 
in the convecting phase and to thermal fluctuations in the 
conducting one, even though in this last case the height 
of the peak is several orders of magnitude lower than in 
the convecting one. Calculations have been performed 
for different system sizes in the patterned region, where 
the peaks are all seen to be centered at approximately the 
same value of k, k ~ fco (in this case fco — 0.63), but their 
heights increase and their widths decrease with increas- 
ing L (see Fig. ^). In fact, it can be seen (Fig. ||) that 
Smax ~ L and AS ~ L , so that the total area under 
the curve S(k) is kept constant with varying L. This re- 
sult agrees with the lack of finite-size effects in J st , since 
the convective heat flux can be seen to be equal to the 
area under the structure function. All these facts lead 
to the conclusion that the structure function approaches 
a delta function in the thermodynamic limit (where the 
rolls are perfectly shaped) as the system size increases. 



where the components of fc have the form 2tt/i/L, where (i 
is an integer and L is the corresponding linear dimension 
of the system. Then it can be seen that the structure 
function is: 



(2.17) 
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FIG. 5. Logarithmic plots of the maximum and width of 
the structure function against system size. Linear fitting gives 
exponents ~ +1 and ~ — 1, respectively. 



III. LINEAR STABILITY ANALYSIS OF THE 
STRUCTURE FUNCTION 

A linear stability an alysis of the homogeneous u = 
solution of model (1.3) permits to estimate the transi- 
tion point as that at which this solution becomes unsta- 
ble. Beyond it the system departs from zero and becomes 
saturated by the nonlinearity, giving rise to the appear- 
ance of a non-zero roll solution. Due to the relevant role 
of the structure function as a means of characterising a 
patterned state, as has been discussed in the previous 
section, we shall begin our theoretical approach to the 
problem with a linear stability analysis of this function. 

We now assume that the homogeneous conducting 
U = solution is affected by small perturbations, so that 
we only need to evaluate the evolution equation of S^, 
defined in (2.17) in the linear regime. This can be done 
by using the Fokker-Planck equation governing the be- 
haviour in time of the probability density of the field (in 
this case in Fourier space). As shown in the Appendix, 
the Fokker-Planck equation for the probability distribu- 
tion of the stochastic process tt^(t) is 



dP 
~dt 



d 



l-i . V 



d d 

, - , T. Gu-V' D V I — Gy + y' P (3.1) 

V ou u ou v 



- E 



where D v is the discrete Fourier transform of -Di-j, and 



F and G are the Fourier transforms of the correspo ndin g 
terms in the discretized version of the SH equation ( [1.3] ) : 



dui 
~~dt 

So that 



= fi(u(t)) + m{t) + 9i (u(t)) &(t) (3.2) 



«-> F M = (r - (Vj + kl) 



g i = U i <r-> Gf, = Up 



(3.3) 



where and V M are a discretized version of the Lapla- 
cian and its Fourier transform, respectively. The time 
evolution equation for is given by: 



dS» 1 /■_ OP 

— = — / LJU U u U-„ —— 

dt V dt 



(3.4) 



When Fokker-Planck equation ( |3.l|) is introduced in 
this expression and integration by parts is performed, 
the following equation is obtained: 

-h— / Du S^ 7] 8_^ y G n _ v < D y f G v + V i P 



v~ v n r p 



\-y J Vu ^2 b-m S/ju, G v -, y < D u , G v+U > P 



-fir] ^fi 



du u 



+ 2eV (3.5) 



Introduction of (3.3) into this expression leads to 
d 



dt 



S, 



'„(*) = 2(r-(v2+fc 2 ) 2 )^(t) + 2 £ + 

S » ^2 D "+ yi ^2D v {u^ v u^ ll+v ) (3.6) 



V 



But, according to definition (2.16), the following relation 
holds: 

J2D V = V D{x = 0) = V D(0) (3.7) 

u 

Thus the equation for the structure function is finally 



6 



p^t) = 2(r + D(0)-(vl + k 2 Q ) 2 )s,(t) 



(3.8) 



Translation of this equation into its continuum version 
leads to the final expression for the evolution equation 
of the structure function of the Swift-Hohcnbcrg model 
in the presence of a space-colored multiplicative external 



^S{k, t) = 2 [r + D(0) - (k 2 - klf\ S(k, t) 



+2e + 2 



1 



(2^ 



D(q)S(k-q,t)dq (3.9) 



The stability analysis of this equation needs some com- 
ments. The last term in Eq. flO| ) is a mode-coupling 
term, which is essentially non-linear and will be discarded 
in our linear analysis. In this case, the presence of a 
multiplicative noise leads to the existence of an effective 
noise-dependent control parameter, r + D(0), so that any 
perturbation of the homogeneous state will always grow 
if the condition r + D(Q) > is obeyed. Hence, lin- 
ear analysis predicts that, under the presence of multi- 
plicative noise, the system can leave the homogeneous 
state in situations for which r < 0. The nonlinearity will 
then stabilise the system in an ordered state. Besides, 
D(0) = D/(4X 2 ) for the particular noise spectrum we 
have chosen, so that the amount of the shift decreases 
with increasing correlation length. This effect can be ob- 
served in the simulation results presented in Fig. |^, where 
a phase diagram of the system is plotted in the (D, — r) 
plane for different values of the correlation length of the 
noise A. It is worth noting that, as long as D ^ 0, the 
critical value of r for the transition from conduction to 
convection is strictly negative. As predicted by our anal- 
ysis, the critical curve in this plane is a straight line whose 
slope decreases as A increases, reducing the noise-favored 
region and thus the shift effect due to the multiplicative 
noise. In the figure, isolated symbols are transition points 
as obtained from simulations, whereas solid lines are the 
corresponding relations coming from the linear analysis 
(r + D/{4\ 2 ) = 0). 




FIG. 6. Phase diagram of the system in the (D, —r) plane 
for different values of the correlation length of the noise. 

In the lin ear stability analysis of Ref. p6| , the last 
term of Eq. ( |3.9|) was considered and was proved to give 
corrections of order D(0) 2 . On the other hand, the rele- 
vance of this term can be seen by studying the s tati onary 
structure function. From evolution equation (|3.9| ), this 
function can be analytically found in the subcritical case 
(r+ < 0) for a white external noise (D(k) = D), the 
result being: 



Sst(k) = 



{ k *-klf+\r eff \ 
where the effective control parameter r e // is 

D 



r eff 



Ax 2 

and the renormalised additive-noise intensity is 

£ 



e = 



1-jD 



(3.10) 



(3.11) 



(3.12) 



with 



7 = -— | r eff | 1/2 arctg \ r eff | 1/2 



1 



(2n) 2 



1 

W) 



dk 



(3.13) 
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continuum space with numerical results coming from a 
discrete-space simulation. Figure ^ shows a comparison 
(for the case D = 0) between these t wo r esults and a (dis- 
crete) numerical integration of Eq. (3.9). The agreement 
between these results and the simulation ones is quite 
good even at large wavenumbers. 

It should be noted that making 7 = in ( |3.10 )-(|3.n 



amounts to ignoring the mode-coupling term in Eq. (f 
Figure shows that the result (3.1C) with 7 = is al- 
most identical to the solution of the full linear analysis, 
in agreement with our previous statement concerning the 
lack of influence of the mode-coupling term in a first- 
order linear stability analysis. 



FIG. 7. Simulation results for the spherically-averaged 
structure function for a conducting state with (circles) and 
without (diamonds) multiplicative noise. Dashed lines are 
the continuum analytical solutions coming from linear anal- 
ysis, and the dotted line corresponds to the solution with no 
mode-coupling term (7 set equal to 0). Values of the param- 
eters are r = -0.2, e = 0.001, Ax = 1.0 and k = 0.63. 

0.006 I . 1 . 1 . 1 



IV. LINEAR STABILITY ANALYSIS OF N-TH 
ORDER MOMENTS 

Due to the stochastic character of the field u (x, t) , the 
linear stability analysis can also be done on the statistical 
moments of u. In this case, differences between the zero- 
dimensional and the spatially-extended cases appear, as 
will be discussed in what follows. 




FIG. 8. Comparison of the spherically-averaged structure 
function obtained by simulation (diamonds), discrete integra- 
tion of the analytical linear evolution equation (dashed line) 
and linear analytical continuum solution (solid line). 

This linear result can be compared satisfactorily with 
(nonlinear) numerical simulations, since in the subcriti- 
cal region nonlinear terms are expected to be negligible. 
Figure [?] shows both kind of results with and without mul- 
tiplicative noise. The ordering role of this external noise 
is observed even in this conducting situation, enhancing 
the maximum value of the structure function at k = k . 
The growing discrepancy between theory and simulation 
at large wavenumbers can be understood from the fact 
that we are comparing an analytical result derived in a 



1. Zero-dimensional case 

Let us consider the following zero-dimensional model, 
known as the Stratonovich model: 



— = ai-i 3 + i^(i) 
<e(i)£(i')> = 2 £>£(*-*') 



(4.1) 



In the absence of noise, this model exhibits a supercrit- 
ical pitchfork bifurcation at a = 0. For negative values 
of a, the stable stationary solution of (4.1) is x = 0. 
In the presence of the multiplicative noise, position of 
the bifurcation point can be given by a linear stability 
analysis of solution x — 0. Since x is now a stochastic 
process, the stability analysis is performed on its statis- 
tical moments, whose time evolution can be found from 
the Fokker-Planck equation governing the evolution of P. 
In the linear regime we have, 



dP 
dt 



d 

■— ax P 
ox 



D^x^xP 
ox ox 



(4.2) 



where P(x, t) is the probability density of the stochastic 
process x(t). This equation leads to 



— (a;™) = n{a + nD) (x n ) 



(4.3) 
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which indicates that the bifurcation point for the n-th 
order moment is located at a c — —nD. Hence the posi- 
tion of the bifurcation point depends on the order of the 
statistical moment which is being analysed, which makes 
this analysis meaningless. We will see shortly that spa- 
tial coupling in extended systems introduces changes to 
this situation, leading to a useful analytical result. 



2. Spatially- extended case 



We now turn our attention back to the SH model with 



multiplicative noise ( |l.3[ ). In order to obtain a possible 
higher-order generalisation of the structure function, let 
us consider the following correlation function: 

G^(x,t) = i J dS '(f[ u(x' + (l-l)x,t)\ 

(4.4) 

which measures the correlation between the values of 
the field at n points equally spaced along a line. Its 
Fourier transform leads to a generalised n-th order struc- 
ture function, 



Hi Hi Hn \i=l / 



(4.5) 



Applying to this quantity the procedure described 
in the previous p arag raphs, making use of the Fokker- 
Planck equation ( |3.l| ), performing integration by parts 
and ignoring internal noise, one finds 

pi n) (t) = n(r + D(0)-(Vl + k^) 2 )s^(t) 



+ 



n(n — 1) 
V 



5>„S£1 (4.6) 



V. CONCLUSION 

Simulations of the Swift-Hohenberg model with a 
space-correlated fluctuating control parameter show that 
the bifurcation point from a homogeneous to a struc- 
tured state is shifted from the standard additive-noise 
case. The amount of the shift increases with increasing 
intensity of the multiplicative noise, whereas space corre- 
lation in the fluctuations shifts back the transition point 
towards the standard situation. A good qualitative and 
quantitative agreement is obtained between this numer- 
ical study and a linear stability analysis of the structure 
function of the system, for the values of the parameters 
used here. The stability analysis has also been performed 
on a set of generalised structure functions which corre- 
spond to higher-order statistical moments, and the same 
results are found in all cases (i.e. for all orders) within 
our assumptions (i.e. at first order), which constitutes a 
qualitative difference to the role of multiplicative noise 
in zero-dimensional systems. 

In conclusion, in this pattern- forming spatially ex- 
tended system, multiplicative noise produces an unam- 
biguous shift of the transition in the ordering direction. 
It is worth saying, however, that according to recent stud- 
ies on other extended systems in the presence of multi- 
plicative noise ]Ti| , one can expect that for large values 
of the intensity of this noise, the system will return to 
the disordered state. 
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Unlike the zero-dimensional case, all dependence on n 
of the linear coefficient in the previous equation is fac- 
torised out. Hence, this results shows that spatial cou- 
pling prevents the linear analysis from giving different 
first-order results for different-order statistical moments 
(part of the contribution of the multiplica tive noise comes 
through the inhomogeneous term in Eq. (|4.6| ) , which rep- 
resents the coupling between spatial modes). There exist, 
of course, other possible methods to study and extract in- 
formation from the behaviour of higher-order moments in 
pattern-forming systems (see, for instance p3| ). 



APPENDIX: FOKKER-PLANCK EQUATION 
FOR A SPATIALLY-EXTENDED PROCESS 
WITH MULTIPLICATIVE NOISE 

In the following we shall assume a multiplicative noise 
with a non-delta correlation only in space: 

( t) t')) = 2 D(x -x") 6(t - t 1 ) (Al) 

Our aim now is to find the Fokkcr-Planck equation 
governing the stochastic process whose evolution is given 
by the general equation 
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du (x, t) 
dt 



= / {u (x , t) , V) + g {u (x, t) , V) £ {x, t) . (A2) 



This can be easily done in Fourier space. In a first step 
the equation is written in a discrete space, where the 
Langevin equation is: 



du 
~dt 



l - = fi{u{t))+ gi (u(t))&(t) 



(A3) 



where the cells have been named with one index inde- 
pendently of the dimension of the discrete space. The 
correlation of the noise in this discrete space is: 



(Si(t)^(t / )) = 2A-i«(*-0 



(A4) 



where Di-j is the discrete version of the continuous func- 
tion D describing the space decay of the correlation. In 
the limit of zero correlation length (white-noise limit) it 
becomes eSij/(Ax) d , where Ax is the spacing of the lat- 
tice, i.e. the cell size. 

Now we define the discrete Fourier transform and an- 
titransform as 

u„ = (Ax) d Y,e-^ x 'u t 

i 

where L is the number of cells per dimension, — j^;P, 
Xi — iAx and the sums are d-fold going from 1 to L. It 
can be seen that the following relations hold: 



i 



The correlation of the noise in Fourier space, £^(t), can 
be seen to be 

< Ut) Ut') > = (LAx) d 2 Dp 5^, v 6(t - t') (A7) 



Here definitions ( A5) and properties ( |A6| ) have been used. 

The equation of evol ution of comes from the Fourier 
transformation of Eq. (A3): 



du 



9i i% 



(A8) 



The last term at the right-hand side can be evaluated 
in terms of the Fourier variables by means of (A5) and 
(P): 



dlLn 

dt ~ U 



(LAx) d ^ ^ 



(A9) 



In the phase space of these Fourier variables we consider 
an ensemble of systems corresponding to a given real- 
ization of the noise and different initial conditions. The 
density of this ensemble p (u, t) must verify a continuity 
Liouville equation: 



dp 

dt 



y — 

^ du„ 



Up p 



(A10) 



where p (u, t) = (S (u(t) — u)) IC , the average being taken 
over initial conditions only. On the other hand, the aver- 
age of this density over the noise is the probability den- 
sity of the stochastic process P (it, t) = (p) (Van Kam- 
pen's lemma [||J ) whose evolution equation is the Fokker- 
Planck equation we are looking for. Performing this noise 
average on (A10) leads thus to 



(All) 



dt ^ duu, (LAx) d ^ du. 



The average in the last term of the second member of this 
equation can be evaluated by means of Novikov's theorem 
p5| , which states that for a gaussian stochastic process 
£ the following relation holds: 



(p («,*)&(*)>= fdt'Y^ium^)) 

Jo 



dp 



e£„(tf) 



(A6) ana ^ usm g (A7) one finds 



(p{u,t)£ v (t)) = (LAx) d D„ 



dp 



(A12) 



(A13) 



And this last average can be calculated in the following 
way: 



dp 



E(MA-'i«iO-«: 



5£- v (t) du n (t) 



y— 

rim.. 



du v \ 6£- v {t) 



5{u{t)-u)) (A14) 



u{t)—u 



The functional derivative in this expression can be cal- 
culated from Eq. (|A9|). The result is: 
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Su v (t) 1 
8£-v(t) ~ (LAx)' 



■9r,+v («(*)) 



(A15) 



so that the Fokker-Planck equation we are looking for is 
finally 



dt 



+ 



Id d 
V^W E ^7 5m- A, 7^ .9,+, P (A16) 



(LAx) 



In the particular case of a noise which is also white in 
space: 



$i,0 



(Ax) d 

and the Fokker-Planck equation is 

---Y — fp 



D„ 



(A17) 



du. 



(LArV ^ L» g »- V ik 9 ^ P (M8) 



E 
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